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Abstract. We compute the two loop MS correction to the Gribov mass gap equation in the 
Landau gauge using the Gribov-Zwanziger Lagrangian with massive quarks included. The com- 
putation involves dilogarithms of complex arguments and reproduces the known gap equation 
when the quark mass tends to zero. 
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1 Introduction. 



The quantum field theory underlying the strong nuclear force is Quantum Chromodynamics 
(QCD). It is an extension of Quantum Electrodynamics (QED) where the gauge fields are 
required to be elements of a non-abelian colour group, SU{3), as opposed to the abelian U{1) 
of electric charge. Whilst this is a simple mathematical generalization the properties of the 
Yang- Mills field theory are significantly different. Clearly QCD is asymptotically free which is 
not unrelated to the fact that the basic fields analogous to electrons correspond to particles 
which are never isolated in nature, called quarks. They are held together in pairs or triplets 
by the quanta of the strong force called gluons. Equally these have never been seen isolated 
in experiments but rather at high energy they are effectively massless asymptotically free fields 
which to all intents and purposes behave as massless fundamental particles. To a degree this 
behaviour is parallel to the properties of the photons and electrons of QED. However, both 
fundamental forces differ in behaviour in the infrared region. For instance, in QCD infrared 
slavery dominates the confinement picture and the gluon propagator does not have the behaviour 
of a massless fundamental particle. One situation where this property can be manifestly seen 
is in Gribov's construction of the gluon propagator at low energy in the Landau gauge, [1] . An 
additional divergence in the structure of QED and QCD emanates from the way one tries to fix 
a (linear) covariant gauge. In QED one can fix the gauge in a global sense. By contrast, Gribov 
pointed out, jl], that in Yang-Mills theory the covariant gauge condition for the Landau gauge 
has an ambiguity. This occurs at zeroes of the Faddeev-Popov operator when different gauge 
configurations satisfy the same gauge fixing condition. In a local region in the neighbourhood of 
the origin of configuration space, where perturbation theory is valid, there is no such ambiguity 
and standard perturbative calculations are perfectly adequate to describe ultraviolet behaviour. 
However, to properly fix the gauge globally the problem of Gribov copies must be taken into 
account in defining the path integral of the theory, [Ij. Gribov achieved this by restricting the 
path integral to the region of configuration space containing the first Gribov region, denoted by 
0. This is defined to be the region containing the origin where the Faddeev-Popov operator, 
M{A) = — d^D^{A), is strictly positive. Consequently, the path integral is cutoff and a natural 
mass parameter, 7, called the Gribov mass emerges, |1J. It is not an independent parameter 
of Yang-Mills but is non-perturbative and satisfies a gap equation. In turn this gap equation 
derives from the restriction of the path integral to $7 by the no pole condition, [1]. In other words 
the average of \/M{A) over is finite. This construction radically alters the infrared properties 
of the theory. For instance, it leads to a gluon propagator which is not fundamental in the sense 
that it has no (real) pole, [1]. Moreover, it is suppressed in the infrared since it vanishes in 
the infrared limit. Further, the gap equation implies that the propagator of the Faddeev-Popov 
ghost is not fundamental but has a dipole behaviour at low momenta which is referred to as 
ghost enhancement. These infrared properties of the constituent fields are believed to be related 
to confinement, [1], and over the years has led to intense interest in studying gluon and ghost 
2-point functions on the lattice and with Dyson Schwinger equation (DSE) methods. 

Another approach was also developed, however, in a series of articles by Zwanziger and 
collaborators, [2l[3llll[5l[6l[7lEl|9l HO], with other relevant contributions in, for instance, 
[111 I12j . In essence the semi-classical approach of Gribov for Landau gauge Yang-Mills was put 
on a firmer footing with the construction of a localized renormalizable Lagrangian, O HI [71 E] • 
The renormalizability being established by various authors, [U [HI [14]. The implementation 
of the horizon condition defining Q in the original approach led to a non-local operator in the 
action which clearly inhibits direct calculations. In [31 [U [Tj H] Zwanziger localized the non- 
locality with a (finite) set of extra fields which defined the horizon condition in an equivalent 
fashion. The beauty of the renormalizability, [8lll31[T4]. aside from allowing for calculations was 
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to demonstrate that none of the known and accepted properties of QCD at high energy were 
changed or upset. For instance, asymptotic freedom remains with the same /3-function. However, 
the advantage of the new formulation was to allow for loop calculations and the extension to the 
next level of computation of the gap equation, gluon suppression and ghost enhancement. This 
was achieved in [Lb} and . In the former the two loop MS gap equation for 7 was established 
when massless quarks are present. This was checked in a non-trivial way by verifying that ghost 
enhancement was satisfied at two loops precisely when 7 obeyed the gap equation. Indeed the 
theory has no meaning as a gauge theory unless 7 does this and hence is not an independent 
parameter of the theory, In the latter article, p^, the one loop gluon suppression was 
verified as well as the exact evaluation of all the one loop 2-point functions of the fields of the 
Gribov-Zwanziger Lagrangian. 

Given this background we come to the main purpose of this article. Clearly in the real 
world quarks are not massless but massive. Therefore, to have a more realistic understanding 
of the Gribov situation it seems appropriate to include massive quarks. As will be evident from 
what is recorded here this is far from a trivial task. First, quarks only appear diagrammatically 
in the gap equation first at two loops. Moreover, this results in Feynman integrals involving 
three scales. Aside from the quark mass itself, the gluon propagator actually has two mass 
scales in the sense of a conventional fundamental propagator. These are ^i^/Ca"/'^ where the 
mass is actually imaginary. (The presence of \/Ca stems from our conventions which follow 
those derived in I16j.) The presence of the imaginary mass further complicates Feynman 
integral evaluation since some of the fundamental functions of one and two loop integrals, such 
as dilogarithms, need to be considered for complex arguments. Therefore, it is the main purpose 
of this article to extend the massless quark two loop MS gap equation of [H] to the massive 
quark case. Moreover, we will discuss the effect it has on the enhancement of the Faddeev- Popov 
ghost. Finally, we note that given recent developments concerning the scaling versus decoupling 
solutions, [Ul [IHl [ini [201 ED 122] , for which there has yet to be a definitive resolution, we note 
that our computations will be the foundation for extensions to the decoupling gap equation. 
This will be required if that solution is eventually established as the correct picture. Moreover, 
this is possible in our approach because the decoupling solution can be accommodated in the 
Gribov-Zwanziger formulation, j23l I24j. Though it will in fact be a more difficult task than the 
current work due to the generation of mass for the localizing Zwanziger ghost fields. 

The paper is organised as follows. Section two is devoted to reviewing the relevant aspects of 
the Gribov-Zwanziger formalism for the massive quark two loop gap equation. The construction 
of the two loop scalar master integrals to the finite part is presented in section three where 
we discuss at length their expression in terms of functions of real variables. This is necessary 
in order to produce a real gap equation rather than a form which has functions of complex 
variables due to the gluon widths. Our main result is provided in section four whilst we draw 
our conclusions in section five. 



2 Formalism. 

In this section we recall the relevant aspects of the basic Gribov-Zwanziger Lagrangian we will 
use to extend the results of . From [31 [U [71 ISj the (bare) Lagrangian is 



^GZ ^ ^QCD ^ ^'^^i'd" {Dy>p^Y^ - u^^^^'d" {D^UJ^ 



ab 

T ( sabc Aau±bc sabc Aau2bc\ dN^'J 



gpbcQU-ae ^jj^^^b ^ecf, ^ 1^ (fabc ^bc _ pbc^af^^bA _ (2.1) 

v2 ^ / 2(7 
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where we use the usual hnear covariant gauge fixing prescription 



= - tG^^G"'^^ - — - c^d^'Dac'' + i^jJ'^m'^ - ruai^'^il^'^ . (2.2) 

Although we will work strictly in the Landau gauge we have included the usual gauge fixing 
parameter a since it is required to derive the gluon propagator. Aside from this it should be 
understood that a is set to zero throughout. Briefly our conventions in (j2.ip and (j2.2p are 



that A'^ is the gluon, is the Faddeev-Popov ghost, is the quark with mass niq and (j)*^ , 
(j)'^^ , u!^ and are the Zwanziger localizing ghosts. The latter pair are anti-commuting like 
the Faddeev-Popov ghosts whereas t/)^'' and 0^'' are commuting. The Lagrangian is expressed 
in d-dimensional spacetime since we will use dimensional regularization throughout to isolate 
the divergence structure of the Feynman graphs where d = 4 — 2e and e is the regularizing 
parameter. The various indices have the ranges l</<A^,l<a< A'^ and 1 < i < Np 
where Nf is the number of quark fiavours and Np and A'^ are the respective dimensions of the 
fundamental and adjoint representations. The various covariant derivatives are 





= 5^c" 











{D,cl,,r = - 9r''Al<l^f (2.3) 



where g is the coupling constant, GJ^^ is the usual gluon field strength and T° are the generators 
of the colour group which has structure functions f'^^'^. We note that we have reverted to the 
conventions of the original form of the Lagrangian, [HIS], in the mixed 2-point sectoiEl. With this 
formulation of the Gribov-Zwanziger Lagrangian we have checked that the results of the massless 
quark gap equation at two loops and Faddeev-Popov ghost enhancement correctly emerge. The 
fields (p^^ and (j)'^ correspond to the localization of the Gribov horizon condition which originally 
was 

^SWa^^^-t-)) = (2.4) 



and now equates to 



r^'{A-^{x)<i>i'{x)) -- ^^^^ 



2 



V2g 

r'^{A^^{x)^'';:{x)) = -^r- (2.5) 

Our conventions are actually crucial to reproducing the correct form of the gluon propagator 
of the original Gribov article, [1]. Using other conventions could lead to, for example, a gluon 
propagator which has a normal mass as well as a tachyonic mass. From (j2.ip and (j2.2p we have 
checked that the propagators of the fields, with momentum p, are 

*In [151116] the Feynman rules of this Lagrangian were used within the computer algebra computations though 
the actual Lagrangian recorded in the articles followed the conventions of [14) . 
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in the Landau gauge where 



[p^ + 771,2] 



p2 



(2.6) 



(2.7) 



is the usual projector. We have retained a non-zero a in inverting the matrix of 2-point functions 
in the quadratic part of the Lagrangian before setting a = to recover the Landau gauge. The 
Feynman rules for the vertices have no convention complications and are straightforward to 
derive from (j2.ip and (|2.2p . Though we note that the explicit cubic interaction of (j2.ip is 
completely passive since it is never present within the Feynman diagrams contributing to any 
Green's function of interest at the two loop level of this article. 



As (|2.ip is renormalizable and incorporates the Gribov properties we now discuss the set- 
up for our computation. The gap equation satisfied by 7 is defined by the no-pole condition 
determining the boundary of Jl, [1]. In the original approach of [1] this equated to evaluating 
the vacuum expectation value of /"^'^j4"^(/)^'' and ensuring it satisfied (12. Sp where the right side 
is a finite object, [3l[ll[7l[8]. However, at this point we note that in all the vacuum expectation 
values one has to take into account the renormalization of the fields and parameters. In this 
respect we note that the anomalous dimensions of all the quantities we require are available 
at three loops in the MS scheme for an arbitrary colour group, [HI [Ml [I5l [16]. Further, at 
four loops the renormalization of 7 is known for the SU{Nc) Lie colour groups, [25]. These 
follow partly through the renormalizability of ()2.ip , [Sj \TE[ [14] , but also because the localizing 
fields and 7 do not undergo independent renormalization in the Landau gauge. Instead all the 
renormalization constants are determined by Slavnov- Taylor identities, [U [T3l I14j . Denoting 
the associated anomalous dimensions of a field or parameter F by 7r(o) where a = /{IQtt^) 
then we record the renormalization constants we require as being encoded in the anomalous 
dimensions 



lA{a) 
77(a) 



[STpNf -13Ca]- + [iOCATFNf + 32CFTFNf-59Cl 



+ 0(a3) 



:CAa + 



AOCATpNf - 95Cj 



48 



[WTpNf - 35Ca] ^ + [280CArpiV/ - 449(7^ + l92CFTFNf] — + 0{a^) 
48 L J 192 



with the /^-function 

Pia) = - 



IIca - ^TFNf 



34 



20 



3 3 

The elementary group Casimirs are defined by 



—Ci - ACfTfNj - —CaTfNj 



(2. 



(2S 



Tr ( T^T^ 



TfS 



ab 



rpCLrj-iQ 



acd fbcd 



cfi , r'^^f 



CaS 



ab 



(2.10) 



Whilst the higher order expressions are available we only provide them at two loops as that is the 
order we compute to here. In (j2.5p we note that the renormalization of all fields and parameters 
present is therefore already fixed and hence after all contributing Feynman diagrams have been 
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computed and assembled the resulting vacuum expectation value is finite. With a massive quark 
present, its mass will be renormalized in principle too. However, as it first appears at two loops, 
scaling it from a bare to a renormalized parameter will not affect the two loop gap calculation 
as the counterterms from the quark mass renormalization constant will only arise first at three 
loops. 

As (12. Sp is the vacuum expectation value of two fields it is easy to determine since essentially 
it is the closure of the legs on the mixed propagator of (j2.6p and integrated over the momentum 
p. Thus for higher loop calculations one simply evaluates the relevant Feynman diagrams which 
are merely vacuum bubbles with various configurations of masses. We devolve to a later section 
the more detailed structure of such two loop massive vacuum bubbles and concentrate in the 
remainder of this section on more general aspects of the two loop gap equation calculation. The 
main ingredients are the generation of the Feynman graphs via the Qgraf package, [26], and 
its conversion into the symbolic manipulation language FORM, [27] . We use Form as it is ideal 
for handling the underlying algebra in an efficient manner. For the gap equation, due to the 
mixed propagators, there are 1 one loop and 17 two loop Feynman diagrams to be determined 
exactly as a function of 7 and niq. As they resolve into the basic structure of two loop vacuum 
bubbles, we note that to make contact with known results we apply elementary partial fractions 
to the common factor in the propagators of the Gribov related fields, such as 



Moreover, within our FORM routines the Feynman rules are substituted automatically and the 
elementary group theory is evaluated making extensive use of the Jacobi identity for the structure 
functions, partly due to the form of the pure (j)'^ propagator. Therefore, all that remains in 
determining the gap equation for massive quarks is the substitution of the explicit forms for 
the master Feynman integrals which the Form routines produce. The next section is devoted 
to this where we concentrate on the intricacies of dealing with vaccum integrals with massive 
quarks and complex gluon masses. 




3 Master integrals. 

There are two main scalar master integrals which arise in the computation. The first is the 
massive one loop vacuum bubble which is virtually trivial in comparison with that we have to 
consider at two loops. Though it does arise in the two loop computation when a line in the basic 
form of a two loop vacuum bubble graph is omitted. Therefore, defining 

where 



(3.2) 



includes the d-dimensional momentum space measure we have exactly 

which is trivial to expand in powers of e. Therefore, we now concentrate on the basic massive 
scalar two loop vacuum bubble which we define as 

h{ml,ml,ml) = f f . (3.4) 

^ A Ji [fc2 + ^2] _ 1^2 + ^2 [/2 + ^2] 
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which is completely symmetric in its arguments and has been studied extensively over the years. 
See, for example, |28l [29l I30j . Its expansion in powers of e, where d = 4 — 2e, is known to 
several orders but for our purposes it suffices to record it to the finite part. For this we follow 
the notation and conventions of [29]. Then we have 



(47r)^/2(x,y,z) 



c 

1 



2 



[L2-6L1 + ^{x,y,z) + c(7 + C(2)) 
+ {y + z — x) ln(y)ln(z) + {z + x — y) ln(2;)ln(y) 
+ {y + X - z)Hy)Hx)] + 0(e) (3.5) 



where we define 



xln (x) + y\n (y) + zln (z) 
X + y + z 



+ + — 2xy — 2xz — 2yz 



1/2 



We also use the same notation as ^29j in defining 



ln(m^ 



In 



m 



(3.6) 



(3.7) 



where /i is the mass scale which enters when using dimensional regularization to ensure the 
coupling constant remains dimensionless in d-dimensions. The key part of this e expansion is 
the function ^{x, y, z) whose explicit form depends on the sign of the combination of masses 
denoted by a^. For our purposes we note that for > then, [29], 



where 



and 



C(x,y,z) = 8a[M{(^,) + Mi<Py)-M{-(^,)] 

M{(j)) = - [''' dein{smh{e)) 
Jo 

c — 2x' 



coth' 



Moreover, we note that 



coth"^(z) 



1 



■In 



2a 

z + 1 
z - 1 



(3.8) 
(3.9) 
(3.10) 

(3.11) 



and the integral defined by the intermediate function M((p) can be written in terms of known 
functions 

M(0) = </.ln(2) - + 1^(2) - Li2(e-^) - Li2(-e-*) (3.12) 



where Li2(2:) is the dilogarithm function, 

Li2(z) = 



ln(l - x] 



dx 



(3.13) 



and ^(2:) is the Riemann zeta function. As an exercise to aid the interested reader it is instructive 
to consider the elementary case 12(0, 0,m^) which is explicitly 



12(0,0, m^) 



kJi k^{k-iy[P+m^] 



(3.14) 
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It can be evaluated directly and then compared with ()3.5p to give 

,2 



9 2 

/2(0,0,m2) = - _ - _^3-21n(m2) 



7 + 3C(2) + 21n^(m2) -61n(m2)l + 0(e) 

(3.15) 

which will be required for checking our expressions in the massless quark limit. 



However, as we are ultimately interested in the massive quark case, we have to consider 
several master integrals. These are 



/2(mq,m^, -i\/CA7^) 



1 



k Jl 



A;2 + 



(3.16) 



and the related integrals 



+ rriq 



{k-lf+m? [l^ + i^/Cl-i'^\ 



k^+m? {k-iy+ml [P-iy/Cl-f'^ 



(3.17) 



We concentrate on the former two as the definition of the latter follows from using elementary 
calculus. For the first we will focus on 



i{WCAl\mlml) = 8a 2M(</)„2) -M(-0,^^2) 



where now 



(3.18) 
(3.19) 



leading to the intermediate variables 



coth 



coth^i 




(3.20) 



Though for practical purposes it is more appropriate to re-express these by applying the loga- 
rithm definition 



e 1 



Cai" + \Ca1^ + ^i^Al'^ml 



"2m2 



These naturally lead to the two functions 

VCI7' + JCAl^^^i^Al'ml 



(3.21) 



M(t, 



^ + im 
2 2 



/CI72 - \ Ca1^ + 4iV^72m2 



ln(2) 
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and 



M(-0, 



Li2 

2 







+ 4^^ 






V 


+ 








+ 


fCAl'^rnl 






+ 4i\ 





C(2) 



+ In 



Ca7' + ;/Ca7^ + 4iV^72m2 



/CI7' - v/cuy + 4V^ 



-In^ 



Li2 



+ JCA7^ + 4iVCl7X' 



- ^CaI^ + 4i VCI7 



- Li2 



^/Cl7' 


- ^CaI^ + 4h 




^/Cl7' 


+ ^CaI^ + 4h 


/CU7^"T'g_ 


yCA7^ 


+ ^iyJCAl'^fn^ 






+ JCaI^ + 4i\ 





ln(2) 



(3.22) 



(3.23) 



where we have used the relationship, [31 



Li2(x) + Li2(-x) = -Li2(x^) 



(3.24) 



Remarkably, this leads to the compact expression 



i{iy/c2'^'^ ,w? ml) 



M\ Cai'^ + 4i VCa7^"^o 



C(2) 



+ Tin' 



Cai" + \Ca1^ + 4iVCU72m2 



+ Li2 



C7^72 - ^CA7^ + 4iVCU72m2_ 



^72 + JCaI^ + ^iVCl-i^ml 



(3.25) 



giving the integral to the finite part 



/2(mg,mg,iA/CU7^) = - ^{^VCai"^ + 2m; 



2e2 
1 fl 



-[-{^WCAr +Qml) - 2m^ln(m^) - WCaYH^^ 



^7' 



1 



iVCA7'ln'(iv/CA7") 



+ 6m2ln(m2) + ^WCaYHWCaI 



- 2iJCA7^ + WCAl^rn^, 



C(2) 



+ Tin' 



/CI7' + VCa7' + 4iV^72m2 
/CI72 - yCUy + 4V^7^ 
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+ Li2 



UiVcli^ + 2ml) (7 + C(2)) 



^(2m2-,yCU7') InK) 



iVC^jMml)HiVC^j^) + 0(e) 



(3.26) 



We have checked that this expression correctly reduces to that for l2{0, 0, i\/ Caj'^) in the hmit 
— > 0. This is not as straightforward as it seems due to the presence of the difogarithm 
function and terms involving ln(r?T,g). Disregarding all terms proportional to and expanding 

CaJ^ + 4:^^yCA7'^rnq in powers of m^, then l2{mg,mg,i-yCA'j'^) reduces to 



him m i^CAl ) 



2e2 



Hi^CAY 



WCaY 



C(2) + 2Li2 
1 

+ 2 



im" 



ln"(VCA7^) - 21n(VOA7^)ln(m^) + ln"(m; 



^^^'(7 + C(2)) + 



In^(m^) 



2 V -V // 2 
iv/cl72IK(m2)E(iv/C^7') + 0(m^;e) 



(3.27) 



where we note that the imaginary dilogarithm vanishes as and the remaining logarithmic 

terms in m^ cancel. By making the analytic continuation — > iy/CA'y'^ in (I3.15p . we see that 
our integral him"^, m"^, i\/CA'y'^) is entirely consistent with ^{0, 0, i\/CA7^) in the limit of zero 
quark mass. 



2M{(t>r, 



M{-(t>_ 



Next we turn to the complex conjugate integral and focus on 

^(-^\/CU7^"l^,mq) 

where now the variables are 

i 



a 



CaI^ - U^CAl^ml 



2mt 



(3.28) 



(3.29) 



leading to 



CaY + \CaY + ^i^Yml 



CaY + ^i^Yml 



CaY + ^i^fClYml 



\ VC^Y + ^CaY + ^i^Al- 
Without reproducing analogous manipulations, we find 



2^2 



(3.30) 



i{-i\fC~Al^ ,ml,ml) 



AidCAY - U^CAYml 
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41n(2) 



+ 



C(2) 



+ Li2 
+ 4E 



Ca7^ + 4iVCl7'"^i 
CaI^ + 4iV^72;;^ - ^^72"" 



where we have used 



CaI^ + ^i^Al'ml + ^72 
Cai" + v/CA7^ + 4iv^72m2 

ln(l - z) 



+ 0(e) (3.35) 



(3.36) 



Again we have checked that the correct massless quark hmit emerges with the direct evaluation 
of the equivalent integral. 

Whilst we have now determined all the master integrals to the finite part in the e expansion, 
the explicit expressions are not in a fully useful format. Given that the ultimate gap equation 
is a real function we need to write the expressions as a real and imaginary part. This is not 
a simple exercise due to the presence of the dilogarithm of a complex argument. However, the 
theory behind such functions is known, |31j . and we summarize what we require for the current 
calculation. Writing the complex variable z in polar form we have the real and imaginary parts. 



Li2(re^^) = Li2(r,0) + i [u;ln(r) + iCl2(2w) + iCl2(2e) - iCl2(2u; + 2Q)\ 



where 



Li2(r,0) 



1 r ln(l-2a;cos6l + x2 



-dx 



and 012(0) is the Clausen function defined by 



Cl2(0) 



In 







2 sin 



The intermediate angle oj is related to the polar variables r and of z by 

r sin0 



tan' 



-1 



1 — r cos / 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



Given these general definitions then to proceed with our simplification to real and imaginary 
parts, we need to write the arguments of the dilogarithms in polar forms. To assist this we recall 
the elementary lemma for a complex variable z = a + i6, where a and h are real. 



1 

72 



Va2"+62 + a ± 



V2 



a . 



(3.41) 
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So, for example, 



(3.42) 



For the dilogarithms if we set 



re = 



^Cai^ + M^Al'ml - ^Al" 



(3.43) 



then 



re 



giving 



(3.44) 



r = 



tan^ 



^(7278 + \^CAl^m\ - CaI^ 



(3.45) 



In what follows wc will always regard r and 9 as taking these values with the associated corre- 
sponding value of LO. Although the dilogarithm is the most involved of the terms which appear 
in the finite parts, similar manipulation is required for several of the logarithm terms. Collecting 
all the pieces together we find the following expression written as real and imaginary parts, 

1 /I 



(3^^/CA7' + 6m^) - 2m^hi(rn^) - WCaiMWCaY) 



€ V2 

m2E'(m^) - ^i^CAl'^\i^CAl' 



2 



X 



-In (^Cl7« + l&CAl^m\ + Cai") 

IE i^Vcl^^ + y 78 + 16CA7^m4 + CaI^ 



+ 2ztan"^ 



CI78 + IQCAl^mt - CaI^ 



V2 ^^7^ +s]yJCll^ + 16CA7^m4 + 



1 . ^, 

2' 



C(2) 



2 ln(2) + - H^Cai'') - + ^ 
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+ Li2(r 
+ iu 



21n(2) + -ln(CA7^)+ln(m: 



- IE (^V2Vc2j^ + \j ^Cl78 + 16CA7^m4 + (7^7^^ 



+ ^Cl2(2a;) + ^Cl2(2^)-^Cl2(2a; + 2e) 



X In ||^/2\/Cl72 + y ^Ci78 + 16CA7^m4 + CaI^ 



+ 2ztan"^ 



CI78 + \^CAl^m\ - CaI^ 



- 21n(2) + liTT - 111(^017') - InK) + ^ 



+ Li2(r 
+ io; 



21n(2) + -ln(CA7^)+ln(m2) 
- (^Cl7« + 16C^7"m4 + Cai') 



In |^V2V^72 + ^ ^Cl78 + 16CA7^m4 + (7^7^^ 



+ ^Cl2(2u;) + ^Cl2(2^)-|ci2(2cc; + 2e) 



^(^ycl7' + 2m2)(7 + C(2)) - V2m^-zv^72)ln'(m^) 



^^/CA7'ln(m^)ln(^VCA7') + 0(e) • 



(3.46) 



Whilst this is not truly of the form a + ih since not all terms have been fully multiplied out 
and there are logarithms with purely imaginary arguments, we prefer to leave it in this more 
compact form since, for instance, it is elementary to implement 



ln(iVCA7 



,2\ TTi 



ITT 



(3.47) 



within our FORM routines. This also takes care of the other elementary complex algebra auto- 
matically. 



4 Two loop gap equation. 

Equipped with the basic master integrals we are now in a position to assemble the two loop 
Gribov gap equation in the MS scheme with massive quarks. This requires the evaluation of 



14 



the seventeen contributing Feynman diagrams which without the power of Form would have 
been virtually impossible. Having already discussed the key aspects of the computation, we 
ultimately find 



aCA 



+ a' 



+ a' 



^TFNfm, 



4lo + 



TT 



r 

2017 _ 11097 ^ 65 
768 2048 256^^ ^ 48 ^ ' 



35 /-,^ 4,\2 1137 2057r2 



+ CATFNf 



21n(2) - ^ + ^ln\ml) - ^ln{mlMCAl') 



+ l^ln(C^/)-ln 



C'X^^ + 16CA7^ml + Ca7^ 



- In 



\/2x/CU7' + \I\/Ch^ + l^CAl^m\ + Cai" 



+ 



+ a\ + mCAl^ml + CaI^ 



V272 ) 



^ - i ln2(2) - \ ln(2)ln(m^) - \ ln(2) ln(C^7') 



+ ^ln(2)ln 
+ ^ln(mg)ln 



V2 V^72 + J ^6-278 + IQCaiV^ + Cai' 



78 + 16CA72m4 + CaI^ 



V2VC27^ + JJcIj^ + 16C^74m4 + Ca7^ 



1 ^2 



32 



In (Ca7' 



+ ■Jln(C^74)ln 



- ih^^ 



V2^/CI72 + V ^^17^ + 16Ca7^K + Ca7^ 



+ Jln(m^)ln 



CIj^ + 16CA7^m4 + Ca7^ 



^Cl^^ + 16CA7^m4 + Ca7' 



X In 



- ih^^ 



^(7278 + l&CAl^ml + Cai' 



V2 VCI7' + V \lcli^ + 16Ca7^K + V 



^(7278 + 16CA7^m4 + Cai' 



TT TT 

+ r+32 
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+ a' 



- C(2) - 21n2(2) - 21n(2)ln(m2) - lii(2)ln(C47^ 



+ 



+ 



+ 



21n(2)ln 
21n(2)hl 



ln(mq)ln 



- \\n\ml) - \Hml)HCAl') 



iln'(C^7') + lln(C^7')ln 



V2 V^7' + V 78 + lQCAl^m\ + Cai' 



2^" 



^2^0^72 + J JCI78 + 16Ca7^< + CaI^ 



+ 



+ 



ln(m^)ln 



hn{CAl')ln 



CI78 + IQCAl^ml + CaI^ 



C2 78 + l6CA7^m4 + CaI^ 



- In 



X In 



2 



C72 78 + l&CAl^m\ + Ca7^ 

yjc\-i^ + 16CA7^m4 + CaI^ + 2a;^ - 2Li2(r, 0) + ttu; + 



TT 



CI78 + 16CA7^m4 



V2 



4 2 



1 



- ^ ln2(2) - - ln(2)ln(m2) - ± ln(2)ln(C^7') 



+ 



iln(2)ln 
+ -ln(2)E 



Ci78 + 16C^74m4 + Cai'^ 



1^2 



In (m 



2\ 



1. 



ln(m^)ln(CA7' 



+ 



^ln(m^)ln 



+ -ln(CA7*)ln 
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+ 7^n(m^)ln 



'c\-f^ + 16Ca7^K + ^^^^ 



X In 



'Ci7« + 16C^7^m4 + C^74 + _ ^Li^Cr, 0) + ^ + ^ 



+ a" 



'Cl7« + WC^^^ml - C^7^ (^^^^) 
- J ln(2) - ^E(m2) - ^HCai') + y^Cl7« + 16C7^7M + ^^7^ 



+ -In 



_ 'Cl 78 + 16C7A7^m4 + C7^74 
+ ^012(20) -^Cl2(2^ + 2w) + ^CL 



X 



2(2c^) 

7rln(2) + + jHCai') - ^IH 



'CI78 + 16CA7^m4 + C7^74 



2 



^VCl7' + y \/ch^ + 16Ca7^K + Ca7^ 
Cl2(20) + Cl2(2e + 2io) - Cl2(2w)]] 



'^278 + 16Ca7^"i4 _ CaI^ 
'C2 78 + 16Ca7^K 



^Cl2(20 + 2^) - ^Cl2(20) - ^Cl2(2^) + ^ ln(2) + ^Rm^) 

'C2 78 + mCA-f^m^g + Ca7^ 



+ — ln(C^7^) - -In 



16 

- 1e 



V2VC^j^ + ^J^Clj^ + 16CA7^ml + CA7^ 



+ 0{a^ 



(4.1) 



This is a real expression and the main result of our article. There are several checks. Whilst 
we have been careful in checking that the four two loop scalar master integrals reduce to the 
correct expressions in the massless quark limit, the overall final gap equation must also satisfy 
the same test. We note that (|4.1|) does do this and for completeness note that one obtains 



Ca 



5 3 , 
In 



17 



+ 



2 / 2017 11097 95 65 / ^^7^ ^ 35 / Ca7^ 

W^^+ 256^(')-48^"l,^j + 128(,^" 1,^^ 



2 




(4.2) 

where = (2V3/9)Cl2(27r/3) which was originahy recorded in [15]. However, there is another 
check on (|4.ip which is to examine the Faddeev-Popov ghost 2-point function in the zero mo- 
mentum hmit. As was noted in [T] there ought to be ghost enhancement which equates to the 
Kugo-Ojima criterion being satisfied, [32]. Formally writing the radiative corrections to the 
Faddeev-Popov ghost 2-point function as u{p^) then ghost enhancement follows if u{0) = — 1 
which is the Kugo-Ojima condition. This was verified at two loops in the massless quark case in 
|16j . Therefore, we have repeated that calculation here with massive quarks and examined the 
zero momentum limit. This involves applying the vacuum bubble expansion to the 31 contribut- 
ing two loop Feynman diagrams. The computation makes use of the master integrals discussed 
in section 3 and we have used the same routines in order to do the Form identifications. The 
outcome is similar to |16j . In other words the Kugo-Ojima criterion is satisfied at two loops 
precisely when the two loop massive quark Gribov gap equation is satisfied. Indeed as empha- 
sised in Zwanziger's articles, the theory has no meaning as a gauge theory unless this occurs. 
Therefore, we are confident that our result (14. ID is correct. 



5 Discussion. 



We conclude with several observations. The inclusion of massive quarks in the Gribov- Zwanziger 
approach has not affected the main properties of the Faddeev-Popov ghost enhancement at two 
loops. Moreover, the one loop verification of gluon suppression of [16] is also unaffected with 
massive quarks. This is because at one loop the diagrams involving massive quarks do not 
arise in that part of the matrix of 2-point functions responsible for the vanishing of the gluon 
propagator in the infrared. It is worth noting that our original expectation was that massive 
quarks would not upset these key properties of the Yang-Mills fields. One of the main outcomes 
of the result (j4.1j) is the very much involved form which is clearly due to the multi-scale nature 
of the underlying Feynman diagram. Whilst we have concentrated on what is now known as 
the scaling solution, [iTl [I8l [El [20l [211 [22], rather than the decoupling solution it does serve 
as an indication of what to expect if one were to study the same problem in the latter case. 
For instance, within the Gribov- Zwanziger context, [23[ I24j. the gluon propagator acquires an 
additional mass scale deriving from the condensation of a mass for the Zwanziger localizing ghost 
fields. Aside from giving three scale two loop integrals for Feynman graphs without quarks it 
will result in four scale two loop integrals for the case we studied in depth here. Clearly that 
would be a difficult computation. 
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